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We study Bose-Hubbard models on tight-binding, non-Bravais lattices, with a filling of one boson per unit 
cell - and thus fractional site filling. We discuss situations where no classical bosonic insulator, which is a 
product state of particles on independent sites, is admitted. Nevertheless we show that it is possible to construct 
a quantum Mott insulator of bosons, if a trivial band insulator of fermions is possible at the same filling. The 
ground state is simply a permanent of exponentially localized Wannier orbitals. Such a Wannier permanent 
wave function is featureless in that it respects all lattice symmetries, and is the unique ground state of a parent 
Hamiltonian that we construct. Motivated by the recent experimental demonstration of a kagome optical lattice 
of bosons, we study this lattice at 1/3 site filling. Previous approaches to this problem have invariably produced 
either broken-symmetry states or topological order. Surprisingly, we demonstrate that a featureless insulator 
is a possible alternative and is the exact ground state of a local Hamiltonian. We briefly comment on the 
experimental relevance of our results to ultracold atoms as well as to 1/3 magnetization plateaus for kagome 
spin models in an applied field. 



Introduction. — Much recent activity in condensed matter 
physics has focused on identifying nontrivial phases of mat- 
ter which are inherently quantum-mechanical and cannot be 
understood by perturbing around a straightforward classical 
limit. Proposed systems where such phases might occur in- 
clude strongly correlated electronic systems, frustrated anti- 
ferromagnets, and insulating phases of bosonic lattice sys- 
tems, the last of which is our focus in this paper. Unlike 
fermions, bosons are precluded from forming noninteracting 
band insulators so all crystalline insulators require interac- 
tions: hence the term 'Mott insulators' [1|. Mott insulators 
realized at integer fillings of bosons per site have a classical 
description deep within the insulating state in terms of a fixed 
integer number of particles per site, which means that to find 
nontrivial states, fractional site filling is desirable. 

Theoretically, it has been proven [2, 3| that insulators at 
fractional filling per unit cell cannot be 'featureless': they 
form either crystals with an enlarged unit cell by breaking lat- 
tice translation symmetry, or exotic phases with topological 
order - i.e. phases with emergent excitations that carry un- 
usual statistics. Experiments on cold atoms in optical lattices 
have extensively explored Mott insulating phases on a variety 
of simple lattices HQ. All these are Bravais lattices, with 
one site per unit cell - so that the site and unit cell fillings are 
identical - ruling out featureless states at fractional site filling. 

However, recently, more complicated optical lattices with 
a basis - such as the honeycomb [8] and kagome [9] struc- 
tures - have been created, and the Mott and superfluid states 
in them have begun to be studied. This naturally leads us to 
consider fractional site filling, but integer unit cell filling. Are 
symmetry-breaking or topological order still the only alterna- 
tives? We will study examples where the answer is no, but yet 
quantum fluctuations of bosons must be significant, even deep 
within a featureless insulating phase. 

As an example of why a featureless Mott phase at frac- 



tional site filling can be counterintuitive, consider the kagome 
lattice at a filling of one boson per unit cell, or 1 /3 site fill- 
ing. Symmetries conflict with the usual caricature of a Mott 
wavefunction: the essentially classical picture of a fixed num- 
ber of bosons tied rigidly to each site. Attempting to draw 
such a classical cartoon on the kagome at the given filling 
leads inevitably to symmetry breaking: for example, a uni- 
form (q = 0) state that distinguishes a single site in each tri- 
angle, or the 'V3 x order which enlarges the unit cell 
(Fig. [TJi. Arbitrarily choosing a unit cell, say an upward- 
facing triangle, and delocalizing each boson across this choice 
of sites gives a more quantum-mechanical insulator, but leads 
to a state which breaks point group symmetries - specifically 
in this 'triangular' state, that of rotation by 180° lfT0l (Fig. [TJ. 
Since the unit cell has three sites, which in the tight-binding 
limit are indivisible, one cannot write a state as a product 
of disjoint 'molecules' while respecting the full sixfold point 
group symmetry. Finally, more sophisticated approaches that 
first implement a duality transformation and condense vortices 
of the dual order parameter also break symmetry at this filling 
ifTTIl . Thus, for the kagome (and similar examples) fractional 
site filling ensures that there is no smooth connection to a triv- 
ial insulator. The usual strategies for constructing insulating 
phases always lead to broken symmetry, and to date no fea- 
tureless insulating state has been found. 

Here, we discuss a general scheme to write down wavefunc- 
tions and parent Hamiltonians for a class of such featureless 
quantum Mott insulators, including the 1/3 filled kagome as 
our primary example. Although we are interested in the prob- 
lem of bosonic Mott insulators, we first study the reference 
problem of fermions at the same filling on the same lattice. 
If a band insulator of fermions exists, for which exponentially 
localized Wannier orbitals (WOs) that respect lattice symme- 
tries can be constructed, then we show that a bosonic Mott in- 
sulating state is also feasible. The Mott insulating wavefunc- 
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FIG. 1. (a) Kagome lattice and notation used in the text. Standard 
methods to construct featureless phases break symmetries, e.g. (b) 
'v3 x v3' state and (c) triangular 'molecular orbital' state. 

tions we construct are permanents lfl2ll of these WOs- hence, 
we dub them Wannier permanent wave functions- analogous 
to fermionic band insulators which are determinants of the 
same orbitals. Furthermore, we demonstrate that these Wan- 
nier permanent wave functions are exact ground states of sym- 
metric, local Hamiltonians of the Bose-Hubbard type. As an 
added bonus, the Bose-Hubbard models we discuss can also 
be viewed as describing XXZ spin systems with fixed magne- 
tization (e.g., in an applied field), for which our wave func- 
tions describe fractional magnetization plateaus which do not 
enlarge the unit cell. On the 1/3 filled kagome lattice, we 
first find a tight-binding model which at this filling describes 
a band insulator, and use this to construct a Wannier perma- 
nent wavefunction for the insulating phase. To the best of our 
knowledge, this is the first proposal of a featureless ground 
state for this problem. 

Identifying a lattice like the kagome is subtle from a sym- 
metry perspective: it entails taking the (experimentally well- 
motivated) tight-binding limit. Absent this restriction, the 
kagome, honeycomb, and triangular lattices are indistinguish- 
able as they all have the same space group. Formally, the 
tight-binding limit identifies a particular representation of the 
symmetry, which contains more information than the space 
group. We can understand this intuitively: the kagome is a 
non-Bravais lattice with three sites per unit cell. Although 
solving a tight-binding model entails a specific choice of in- 
equivalent sites to form the crystal basis (which breaks lattice 
symmetry), the final many-body band insulating wavefunction 
for fermions is independent of this choice and respects sym- 
metries. In contrast, choosing orthogonal real-space orbitals 
to build a bosonic permanent is more sensitive: e.g. up- and 
down-triangle permanents are distinct, and both break 180° 
rotation symmetry. Thus, restrictions placed by tight-binding 
and the symmetrization of many-body bosonic wavefunctions 
makes constructing featureless Bose insulators challenging, in 
contrast to the more usual case of fermionic band insulators. 

The wavefunctions we study are in a sense bosonic analogs 



of the Affleck-Kennedy-Lieb-Tasaki (AKLT) states H2J [H 
of spin systems, which are quantum-mechanical paramagnets. 
As noted in lfT5ll . the AKLT insulators are 'fragile' Mott in- 
sulators of electrons, which cannot be adiabatically connected 
to band insulators without breaking a crystalline point-group 
symmetry. Similarly, our wavefunctions cannot be deformed 
into disjoint product states without breaking translational or 
point-group symmetry. However, since we discuss bosons, 
in contrast to |[T5ll we cannot consider a noninteracting limit 
whatever the symmetry. We note previous work lfl6l that 
used power-law localized WOs to study Bose condensation 
and charge-density-wave order in the kagome flat band is un- 
related to the featureless insulators we discuss. 
Fermionic band insulators as bosonic Mott insulators. — 
Given a set of symmetric, exponentially localized WOs gn(i), 
with R a Bravais lattice vector and i a lattice site, the Wannier 
permanent wavefunction for a bosonic Mott insulator takes the 

form \% w ) = IIr^r! )' where w r = X«5rW & I adds a 
boson to the WO at R. \^w) is the exact ground state of a 
parent Hamiltonian given by an interaction between densities 
projected into the WO basis. The correlation function can be 
expressed as a projector over the occupied band, 

(b\bjh w = 9k(i)9nti) = »)■ C 1 ) 

R 

and decays exponentially, (b\bj) ~ e J I , with £ related 
to the localization length of the WOs. Higher-body correlators 
can be similarly computed, and also decay exponentially. 

WOs are constructed by integrating Bloch wavefunctions 
over the Brillouin zone (BZ), with an arbitrary phase choice 
iySk at each wavevector k |17|. For fermions, the fact that 
the different WOs are related to each other and to the single- 
particle Bloch eigenstates by unitary transformations means 
that 'filling' any such complete set of single-particle states in 
the lowest band and then antisymmetrizing yields the same 
insulating wavefunction, as the resulting Slater determinants 
are all unitarily invariant. However, symmetrized permanents 
such as | ^w) are n ot unitarily equivalent for different choices 
of ffR.(i), making the choice of single-particle states to 'fill' 
crucial: permanents for featureless Bose insulators must be 
built from localized and symmetric single-particle states. 

To preserve symmetry, we must be able to fix the choice of 
localized WO unambiguously. For a simple band - one that 
touches no other bands - respecting time reversal and inver- 
sion symmetry lfl8l fixes ip^ up to an overall sign, which is 
fixed by requiring continuity of (p-^ across the BZ, a necessary 
condition for exponential localization [19-22]. Whenever a 
fermionic band insulator exists at a given filling, we can find 
such a set of WOs giving a symmetric insulating ground state 
and corresponding local parent Hamiltonian- a significant and 
hitherto unnoticed aspect of bosonic Mott physics. 
Kagome Lattice. — To use this approach for the kagome lattice 
at 1/3 filling we must find a fermionic band insulator (sepa- 
rated lowest band) on the kagome. We consider the general 
case Hq [{tij}] = — J^i j Ujblbj on the kagome lattice, where 
the Uj are assumed to respect all the symmetries of the lattice. 
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Hq has three bands, and with only nearest-neighbor hopping 
(ti) the two lower bands form gapless Dirac points (Fig.[2]i. 
Mass terms that can gap out the Dirac points break inversion 
or three-fold lattice symmetry ('triangular' and '-y/3 x 
states, Fig. [TJ. However, as second- and third-neighbor hop- 
pings (*2 , £3) are increased from zero, the bands become de- 
generate at K, K' when ^2=^3 = \t\\ thereafter the lowest 
band is simple with the twofold If -point degeneracy trans- 
ferred to the upper two bands (Fig. [2] For a symmetry analysis 
see E3l ). This procedure preserves lattice symmetry, provid- 
ing symmetric WOs for the lowest band. 
Wannier Permanents. — For clarity, we focus initially on a 
fine-tuned point l24l l25l . where t\ = t% = t$ = t/6, and 
the onsite potential ta = t/3; our results persist away from 
this point as long as the lowest band remains simple. We may 
rewrite Hq as Ho = —t J2n ^r^R by defining operators that 
'smear' a boson over a single hexagon (Fig. [ljdepicts d R for 
R= a 2 ) 
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Here, R = mai + na2 lies on the triangular Bravais lat- 
tice, ai, a 2 and a3 = — (ai + a 2 ) are shown in Fig. [I] 
and a = 1, 2, 3 labels inequivalent sites within a unit cell, 
so that a kagome site i = (R, a). Since 7^ 
<5r,rs the d R operators are not canonical bosons and we can- 
not form a Fock space out of eigenstates of d R d R . How- 
ever in momentum space, H = —tf BZ ,^ ? 2 tiqrfq, where 
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N s 1/2 E R ^' q '% with N s the number of sites. 

By a straightforward computation, [d q , d^,] — e(q)<5q !q ' with 

e(q) = 1 + | cos (q ' a i)- Since everywhere in the 

BZ e(q) 7^ 0, its inverse is nonsingular and we may define 
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[e(q)]" 1/2 4> sothati/o = -tf B 

7t 



d 2 q 



:(q)4 d q- 



now 



BZ(2tt) 2 ' 

in terms of canonical bosons [d q , dL] = <5 q , q <. Hq has a sin- 
gle band of dispersing states created by rf q with energy — ie(q) 
that is gapped away from E = 0. In writing Ho in this form, 
we have chosen a specific superposition of states in a unit cell, 
and therefore a specific Bloch band. The two orthogonal states 
in the remaining Bloch bands form a pair of zero-energy flat 
bands (Fig. [2]). H respects all the kagome lattice symmetries 
and its lowest band is manifestly simple, as discussed above. 
A single boson is placed in the WO at R by the operator 

d 2 q 
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Here, the function 

#r(R', a) = 



BZ 



d 2 q l + e lqa ° 

(27T) 2 V^Mi)' 



iq-(R'-R) 



(4) 



characterizes the spatial structure of the WOs (Fig. [2]). It is 
easily verified from Q that the % are canonical boson oper- 
ators, i.e. [wr, w ri ] = <5r,r'. Filling one boson in each WO 
gives the Wannier permanent wavefunction 1 26 1 . 




FIG. 2. (a) Evolution of tight-binding band structure of Ho; for 
nearest-neighbor hopping (top) all bands touch, while at the special 
point ti = = tz (bottom) the lowest band is isolated and is es- 
pecially tractable, (b) WO plotted on the lattice. Size of the circles 
give |<?H,(i)|; blue (red) depict gn(i) > (< 0). Exponential decay 
of the boson correlator {&£&o) is also shown below the WO. 



Finally, each individual WO respects the point-group sym- 
metry, as we can explicitly verify. Translations only map one 
WO into another and therefore a product over all the WOs 
is invariant under the entire space group of the lattice, as is 
an arbitrary sum of such products, which includes a perma- 
nent of the form l^w)- Thus, we conclude that \^w) rep- 
resents a featureless Mott insulator at filling 1/3 [27 1. If we 
vary the parameters away from the special point, the WOs 
will change. Nevertheless, as long as the WOs can be chosen 
to be symmetric and exponentially localized, the correspond- 
ing Wannier permanent wave function remains featureless and 
has exponentially decaying correlations. We thus have a fam- 
ily of Wannier permanents, each characterized by the hopping 
Hamiltonian H . 

Parent Hamiltonian. — We construct a Hamiltonian for which 
the Wannier permanent wave function \^w) is an exact 
ground state. We do this by first constructing a new noninter- 
acting Hamiltonian, Hq, which flattens the dispersion of the 
lowest band while leaving its wavefunctions unaffected. Now, 
the WOs are actually single-particle eigenstates. We tune the 
chemical potential /1 to be just above the lowest band. This 
leads to a massive ground state degeneracy, which is lifted by 
adding appropriate interaction terms that ensure single occu- 
pancy of the WOs, thus selecting \^w) as the ground state. 

We begin with the tight-binding problem Ho, and from it 
obtain the operators w R . Using these, we construct a differ- 



H'o can be ex- 



ent single-particle Hamiltonian H' Q = — V Y^o '< 
where N is the total boson number. Using (p 1, 
pressed as a symmetric hopping model on the kagome lattice, 
whose locality follows from localization of the WOs. Ma- 
trix elements of H decay exponentially (identically to (b\bo), 
Fig.[2]i, and hopping between sites on non-adjacent hexagons 
is < 3% of that between nearest neighbors. It is readily veri- 
fied that H'o may be rewritten as H' = — ^ R [(jt( + V)h R + 
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m(^r +"r )] where n R = w^ior represents the occupation 
of the WO of the lowest band at Bravais lattice site R, and 
n R ' 2 are occupations of the WOs of the two upper bands. 

For — V < [i < in Hq only the n R are nonzero. There 
are many ways to fill WOs of the lowest band with bosons, 
leading to a huge degeneracy of many-body states. We lift 
this by adding to H' an interaction that penalizes multiple oc- 
cupancy of a WO. This is accomplished by a term H{ nl = 
If Er"r("r — !)■ I n terms of WOs, this is the familiar 
Hubbard interaction stabilizing a Mott phase with fixed occu- 
pancy of each WO, but corresponds to a more intricate inter- 
action in terms of the site operators of the underlying kagome 
lattice [23]. Like H' , _ff; nt is local in the sense that interactions 
between sites decay exponentially with their separation. It is 
easy to show that \^>w) is the ground state of Hy/ = Hq+H^ 
for —V < n < min(0, U - V) ||28]]. It is also clear that all 
excitations above the ground state are gapped. It is evident 
that \^>w) is the unique state satisfying these properties on 
the torus, so it is not topologically ordered as it has no ground 
state degeneracy. In summary, we have shown \^w) is the 
unique gapped and symmetric ground state of 

HH, = -y^n R + ^n R (n R -l)- M iV (5) 

R R 

for — V < p < min(0, U — V). An equivalent Hamiltonian 
can be obtained by first flattening the lowest band and then 
projecting a Hubbard interaction into it [ 29—3 Q . Hw is intri- 
cate and presently challenging to realize experimentally (see 
below). Numerical study of simpler proximate Hamiltonians 
that may yield featureless phases is left to future work. 
Experimental realization. — A natural experimental setting for 
Bose-Hubbard physics is in ultracold atomic gases in optical 
lattice potentials [4]. The kagome geometry was achieved re- 
cently in an experimental setup using an optical superlattice, 
and was characterized by studying properties of an atomic su- 
perfiuid with such a lattice \9\. Mott insulating states within 
this optical lattice have been recently observed 11321 . 

Methods suitable for identifying an insulating state with 
fractional filling have been demonstrated for atoms in simple 
Bravais lattices. For example, Mott insulating states at various 
fillings can be identified by observing plateaus in the density 
of a lattice-trapped gas within an inhomogeneous potential 
EE]]. The underlying geometry of the insulating state can be 
probed via momentum-resolved correlations in the measured 
atomic density [5|, a method that can be employed to check 
against broken-symmetry states. 

However, the fine tuning of next- and next-next-nearest 
neighbor hopping in such lattices may be difficult to achieve. 
One approach may be to utilize Feshbach resonances to in- 
crease the interaction strength so that a critical ratio of interac- 
tion to kinetic energies can be reached already in a very shal- 
low optical lattice, where substantial higher-order tunneling 
may still occur. Alternately, returning to the analogy between 
bosonic tunneling and spin models, longer-range interactions 
can be achieved and tuned in lattice-trapped gases of polar 
molecules [33] or Rydberg atoms. Although \*&w) is n °t sign 



positive, a related sign-positive wavefunction - potentially 
with better variational energies - exists in the same phase, al- 
beit without a corresponding parent Hamiltonian [34|. 

We note that there are candidate materials where it is 
believed that the spin physics is captured by the isotropic 
kagome lattice Hamiltonian. By exploiting the connection 
between the Bose-Hubbard model and a spin system in an 
applied field, the wavefunctions we consider here are natural 
candidates for featureless fractional magnetization plateaus in 
such materials - particularly since their nontrivial sign struc- 
ture suggests that they will give good variational energies [35). 
Concluding Remarks. — Although we have focused on a spe- 
cific example of a non-Bravais Bose-Hubbard model, the 
kagome lattice at | site filling, the path to generalizing our re- 
sults is clear. For a lattice with a q-site unit cell, it is possible 
to construct a Wannier permanent at filling 1/q if we can con- 
struct symmetric, exponentially localized WOs for the lowest 
band. The first step is to find a tight-binding model whose 
lowest band is simple, which often requires going beyond 
nearest-neighbor hopping, as the kagome teaches us ll36ll In 
general, whether such a model exists depends on symmetries, 
as discussed above. For example, on the 1/2-filled honey- 
comb lattice, the two bands that touch in the nearest-neighbor 
limit form a two-dimensional irreducible representation; as 
there is no third trivial band to which the degeneracy may be 
transferred, the touching cannot be removed without symme- 
try breaking [23|; and, on non-symmorphic lattices featureless 
phases are impossible except at special integer fillings l37l . 

We note in closing that the kagome lattice WO is maxi- 
mal on the hexagon of sites centered on Bravais lattice site R 
and decays rapidly away from it (Fig. [2]) This suggests that 
the wavefunction = IIr^rIO) obtained by filling 'trun- 
cated' orbitals restricted to a hexagon remains in the same 
phase as ■ Verifying this conjecture requires a numerical 
evaluation of correlations in the state I*?), which is a signifi- 
cant problem in its own right. One can by analogy construct a 
featureless insulating wavefunction on the honeycomb lattice, 
even though the Wannier construction fails E3l l34l . 
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SUPPLEMENTARY MATERIAL 
Symmetry Analysis of Band Touchings 

As our approach relies on the ability to construct 
exponentially-localized Wannier orbitals, it will be important 
to determine when it is possible to deform a tight-binding 
model without breaking lattice symmetry in such a manner 
as to produce a simple band at the requisite filling - or at the 
very least, to know in advance when such an effort will prove 
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futile. 

In general, there are two ways to produce a simple band: 
either split degeneracies or just move energy levels around so 
that the singly degenerate states are all in one band. In this 
appendix we seek to understand when it is impossible to pro- 
duce a simple band while preserving lattice symmetries as this 
will serve as an important restriction on the approach we take 
in this paper. Consider a general tight-binding model on a lat- 
tice with a specified space group. This is comprised of a point 
group G, combined with translations through vectors in the 
Bravais lattice. If the lattice is such that a single center can 
be chosen about which all the symmetry operations in G may 
be defined, then the space group is termed symmorphic; if at 
least two different centers are required, then the space group is 
non- symmorphic. Non-symmorphicity requires that the lattice 
have screws (a rotation about an axis combined by a transla- 
tion along a vector parallel to the axis and not in the Bravais 
lattice) and/or glides (a reflection through a plane combined 
by a translation through a vector in the plane and not in the 
Bravais lattice), i.e. point-group operations coupled with non- 
primitive translations 1 38 1 . 

We now focus on a tight-binding Hamiltonian on the lat- 
tice, and consider the Bloch hamiltonian h(q) at crystal mo- 
mentum q. For arbitrary q, no operation in G leaves q in- 
variant and therefore is not a symmetry of h(q). However, 
at high-symmetry points in the Brillouin zone, a subgroup of 
operations in G, known as the little group at q and denoted 
G q , leave q invariant or translate it by a reciprocal lattice vec- 
tor, so that h(q) is left invariant, i.e. for g £ G q , we have 
g/i(q)<7 _1 = hq. The energy levels at high-symmetry points 
can be chosen to be simultaneous eigenstates of both the en- 
ergy and the generators of G q ; the levels can be classified by 
the irreducible representations of G q that they transform un- 
der, and protection, or lack thereof, of band touchings at high 
symmetry points may thus be understood in terms of these. It 
is to this understanding that we now turn. 

For a symmorphic space group, all representations of G q 
are regular (i.e. non-projective.) In this case, constructing a 
simple band consistent with the lattice symmetry is impossi- 
ble if at least one high-symmetry point q there are no non- 
degenerate levels and every degenerate set transforms under 
an irreducible representation of G q with dimension > 2. As 
this statement is somewhat opaque, we illustrate through an 
example: we will contrast the behavior of the kagome and 
honeycomb lattices. Both lattices have the symmorphic space 
group PQram with a triangular Bravais lattice and point group 
D6h; choosing the center of symmetry to be the center of a 
hexagon in each case, we can choose the point group gener- 
ators to be a mirror reflection about the horizontal axis and 
7r/3 rotation as shown in the figure. The nearest-neighbor 
tight-binding Hamiltonian for both kagome and honeycomb 
yields a spectrum with Dirac points at the K, K' points of 
the Brillouin zone where two dispersing bands touch; kagome 
has in addition a dispersionless flat band at high energy which 
touches the upper dispersing band at the V point. To deter- 
mine whether the Dirac points can be removed without break- 



ing the symmetry, we need to classify the irreducible repre- 
sentations (irreps) of the little group Gk of the i^-point (or 
equivalently, of the K' point). We shall discuss the case of 
the K point; our arguments apply mutatis mutandis to the K' 
point as well. Since a rotation by tt/3 maps K to K 1 , it is clear 
that i?Tr/3 Gk\ only rotations by 2tt/3 leave K invariant, 
so only the generator of threefold symmetry -R27173 = 
is in the little group. The other little group generator is the 
mirror reflection a 2 - which also leaves K invariant - and to- 
gether these generate Gk — D^ . For the kagome lattice, 
it is evident that the rotation cyclically permutes the sublat- 
tice indices A, B, and G (corresponding sites denoted by red, 
blue and green dots respectively in all the figures.) The mir- 
ror reflection interchanges B and G while leaving A invariant. 
Since both these operations may be written (as is evident from 
our choice of unit cell labeling) as a mapping R — > gR, ac- 
companied by a permutation of elements in the unit cell we 
may represent the generators of Gk by the matrices 
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Since both these leave the vector (1, 1, 1) invariant, it is 
clear that the representation they generate is reducible (which 
we could have deduced from the fact that D^h has no 3- 
dimensional irreps) and decomposes into the direct sum of 
a one-dimensional and a two-dimensional irrep. We can 
show that eigenstates from the pair of bands that touch at 
K in the nearest-neighbor limit form a two-dimensional ir- 
rep of Gk, while the flat band eigenstates must form a one- 
dimensional irrep. As result, it is possible to construct a sim- 
ple band on the kagome lattice as we demonstrate explicitly in 
the paper. To do this, we add higher-neighbor hopping terms 
and tune their strength until a threefold degeneracy occurs at 
K; once this degeneracy point is crossed, the upper two bands 
form the two-dimensional representation and are degenerate 
at K, while the lower band transforms under the trivial repre- 
sentation at the K point. The band touching in the honeycomb 
lattice on the other hand is protected by symmetry: since there 
are no non degenerate bands at K that form a trivial represen- 
tation of Gk, there is no way to remove the band touching by 
a procedure analogous to the kagome example. 

Turning next to the honeycomb lattice, we note that the mir- 
ror reflection about the axis interchanges A and B and maps 
R — > 02R. However, the 2tt/3 rotation necessarily mixes 
sites from different unit cells, unlike in the kagome case: we 
can show that this operation takes R — > i? 2 7r/3R but also 
maps (in an obvious notation) A(R) — » A(R + ai + a 2 ), 
■B(R) — > B(R — ai — a 2 ). In Fourier space and at the K- 
point, it is clear that the representation matrices corresponding 
to these generators take the form 

«2*/3= e - i27r /3 J>"2- ^ J (7) 

It is readily verified that these form a two-dimensional irrep, 
and thus the Dirac point degeneracy is protected by symmetry 



FIG. 3. Kagome lattice. The center of symmetry is shown, and 
our choice of sites in a unit cell is shaded, and the basis vectors of 
the triangular Bravais lattice are shown. The little group generators 
referred to in the text are also depicted: (i) rotations Ri^i?, by 120 
degrees about the center; and (ii) reflections o 2 about the mirror line. 




FIG. 4. Honeycomb lattice. We use the same conventions and 
notations as Fig. [3] 

in the honeycomb case. 

A three-dimensional example of a symmorphic space group 
with a protected band touching is furnished by the cubic per- 
ovskite structure, which is a lattice of corner-sharing octa- 
hedra; alternatively we may obtain this by placing sites at 
the bond mid-points of a simple cubic lattice. We may then 
verify that the energy bands at R = (tt, 7T, 7t) transform un- 
der a three-dimensional irrep of the little group, as follows. 
Rather than construct the entire little group Gr, we will con- 
tent ourselves merely with constructing a pair of elements 
and use them to show that the three-dimensional represen- 
tation is irreducible. The first of these is the rotation i?27r/3 
about the (1, 1, 1) direction (see Fig. [5]). This simply sends 
R — > i? 27r /3R, an d cyclically permutes A,B, and C; the sec- 
ond is the rotation by R^/2 about the z-axis; this maps R — > 
i? 7r / 2 R- an d leaves C invariant, while taking A(R) — s- B(R), 



FIG. 5. Cubic perovskite structure viewed as a lattice of corner 
sharing octahedra. We show only the symmetries required to prove 
that Gr^^n) has a three-dimensional irrep. 



■B(R) — > A(R — ai). The corresponding representation ma- 
trices are (at the i?-point) 





1 








1 

I— 1 




• 





: 


, R-it/2 — 


I 





• 

















which share no common eigenvectors. Any representation of 
Gr must contain these two symmetry operations and thus can- 
not be reducible, and therefore the threefold degeneracy at R 
is protected by symmetry. 

For a non-symmorphic space group, representations of G q 
for q on the surface of the Brillouin zone are in general pro- 
jective. As a result, if G q contains commuting elements that 
fail to commute in the projective representation, symmetry 
requires that every level at q be degenerate. For instance, 
we may show that if G q contains inversion plus a nonsym- 
morphic generator, then every energy level is at least twofold 
degenerate at q and therefore a simple band cannot be con- 
structed without breaking lattice symmetry. An example of 
this is the three-dimensional space group Fd3m, which char- 
acterizes the diamond and pyrochlore lattices: we can show 
that the little group Gx of the X-point contains inversion and 
a four-fold screw rotation, which commute. As the screw is 
nonsymmorphic, the representations of Gx are projective and 
therefore every band at X is (at least) twofold degenerate, pre- 
cluding the existence of any simple bands. For a lucid expla- 
nation of such symmetry-required band stickings, we refer the 
reader to ||39pT1l . 
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Parent Hamiltonian 



In the main text, we introduced a parent Hamiltonian for 
the Wannier permanent. Here we clarify a few of its essential 
properties. 

The first subtlety is that we rewrote the noninteracting por- 
tion of the Hamiltonian, H' a = —V J2r w r ?«r— /i^V in terms 
of Wannier orbitals of the upper bands, H' — — X)r[(m + 
U)n R + ^("-r + n^)]. We observe that the terms cor- 
responding to occupying WOs of the upper bands are con- 
tributed only by the chemical potential term, and are essential 
since the chemical potential couples to the total particle num- 
ber N, which is impossible to expand purely in terms of the 
lowest band. It is easily verified that [^ R ,"- R 1,2 ] = 0, and also 

["r >"r'] = ["r »^r'' 2 ] = °- Thus > for -V < n <0,it is 
possible to set the occupation of the upper bands to zero, since 
the remaining terms in the Hamiltonian Hw = H'a + H[ nt do 
not mix the lowest band with the upper two bands. Also, it 
should be clear that the symmetry and localization properties 
of the upper two bands are irrelevant since they enter H' in a 
manner that manifestly preserves the symmetry and is local, 
since all terms proportional to /i simply combine to give the 
total particle number N — YliK^i- 

A second point is that the final Hamiltonian Hw contains 
terms that look 'onsite' on the triangular Bravais lattice, but 
in fact are quite intricate on the kagome. The explicit form 
of Hw can be seen by rewriting the Wannier orbital boson 
operators w R in terms of the original kagome lattice boson 
operators b\. From Eq.(3) of the main text, 

< = E 3*( R '> a ) 6 W = E 9^)b\ (9) 

R',c* i 

where we have reintroduced the compact notation i = (R/, a) 
for sites on the kagome lattice. Using this, we obtain 

ij V ' i 

ij kl 

where the interaction 'form factors' 

R 

/S = E^(%R(i)SR(%R(0 (ID 

R 

are computed from the Wannier orbital wavefunction. Owing 
to the exponential localization of g R , they decay exponentially 
with separation between any pair of sites. 

In discussion above we specialized to the kagome, but 
it should be evident that with appropriate changes can be 
adapted to any simple band. 
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